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Abstract 

Many exactly solvable models are based on Lie algebras. The pair- 
ing interaction is important in nuclear physics and its exact solution 
for identical particles in non-degenerate single-particle levels was first 
given by Richardson in 1963. His solution and its generalization to 
Richardson-Gaudin quasi-exactly solvable models have attracted the 
attention of many contemporary researchers and resulted in the exact 
solution of the isovector pn-pairing within the so (5) RG- model and 
the equal strength spin-isospin pn-pairing within the so (8) RG- model. 
Basic properties of the RG-models are summarized and possible ap- 
plications to nuclear physics are emphasized. 
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1 Introduction 



Symmetry is one of the most important paradigms in modern physics. Any 
Lie group and its algebra have a naturally defined action on a product of 
spaces (representations). Thus, they are very suitable for a multi-particle 
system with an underling symmetry. Usually, this means that the relevant 
operators are well defined for one particle as well as for any number of par- 
ticles. This allows one to find exact solutions to a problem, with a given 
underlining symmetry, by referring to the relevant representation theory of 
the symmetry group at place. As a result many exactly solvable models are 
build using Lie algebra representation theory. A well know examples are 
the so (3) and su(2) rotational symmetry, the Elliott's u(3) symmetry model 
[U El El [7] , the Wigner's su(4) spin-isospin symmetry [17], and many more 
that play a major role in nuclear physics. For example, the 50 (8) and sp(6) 
Ginnocchio models, the Fermion Dynamical Symmetry Models (FDSM), and 
the three dynamical symmetries of the Interacting Boson Model (IBM). 

A nuclear many-body system near equilibrium can be viewed as subject 
to a mean field Harmonic Oscillator (HO) potential: 
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Since the symmetry group of the 3-dimensional HO is u(3) [4] one can easily 
see its relevance in the description of nuclei. It is well know that one can 
understand the magic numbers and the shell structure of nuclei within the 
3-dimensional HO approximation [llj. Using the HO single-particle states 
one can write a general Hamiltonian with one- and two-body terms: 

H = ^Eidfcii + - ^ V ijik iafaja k ai. 

Here, and are fermion annihilation and creation operators, £j single- 
particle energies, and Vij t ki = (ij\V\kl) two-body interaction matrix elements. 

The simplest extension of the HO Hamiltonian is to add quadruple in- 
teraction terms Q ■ Q and/or spin-orbit interaction L.S. This has been well 
studied by Elliott and his collaborators [U \5\ El [7J. 

Another important interaction is the pairing interaction in nuclei: 

H P = J2 2^71* ~gJ2 a ti a h a t,i a l,r (!) 
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Here rij is the number operator for pairs. The exact solution of the pairing 
interaction between identical particles in non-degenerate single particle lev- 
els was first given by Richardson [HI [15] . His solution and its generalization 
to Richardson- Gaudin quasi-exactly solvable models (RG-models) have at- 
tracted the attention of various contemporary researchers - resulting in the 
exact solution of the isovector proton-neutron pairing in nuclei within the 
50 (5) RG- model [3] -' to be discussed in Sec. 14.11 and the equal strength 
spin-isospin proton- neutron pairing within the so (8) RG- model [12] to be 
summarized in in Sec. 14.21 Basic properties of the integrable RG-models are 
summarized in Sec. [3] and their possible applications to variety of nuclear 
physics models are emphasized in Sec. El In the next section we briefly 
discuss few dynamical symmetry models of importance to nuclear systems. 



2 Some Exactly Solvable Nuclear Models 

A quantum system has a dynamical symmetry if the Hamiltonian can be 
expressed as a function of the Casimir operators of a subgroup chain. A 
typical example of a rank two dynamical symmetry is the Elliott's su(3) 
model that is used in the description of deformed nuclei [U El El H] • 

H def = eN + X Q-Q. (2) 

This Hamiltonian can be rewritten as a linear combinations of the Casimir 
operator of the su(3) Lie algebra = (Q ■ Q + 3L 2 )/4 involving the 

quadrupole-quadrupole interaction Q ■ Q and the Casimir operator L 2 of the 
50 (3) subgroup of the angular momentum. Thus we have a group chain re- 
duction so (3) C su(3) that provides exact solution to our initial Hamiltonian 

#3o(3)csu(3) = eN + ^L 2 + aCf z) . 

Here N counts the number of particles with energy e of a particular HO shell 
and L 2 lifts the /-degeneracy of a particular harmonic oscillator shell. 

A common feature of the dynamical symmetry nuclear models, is that 
they are all defined for a degenerate single particle levels. Since single particle 
energy splitting breaks the dynamical symmetry, it is usually expected that 
this will prevent the model to be exactly solvable. For example, spin-orbit 
interaction I ■ s lifts the total angular momentum degeneracy j = 1 + 1/2 
and j = (/ + !) — 1/2 and destroys the su(3) symmetry [ID] . Although, the 
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single particle energy splitting breaks the dynamical symmetry it may still 
preserve the exact solvability. The pairing model with non-degenerate single 
particle levels, whose exact solution has been found by Richardson, represents 
a unique example of an exactly solvable model with these characteristics 
[HI [15]. The model is exactly-solvable due to the special extension of the 
relevant dynamical symmetry algebra to a spectral Gaudin algebra. 



3 Spectral Lie Algebras 

A Gaudin algebra Q (g) is an infinite dimensional extension of a Lie algebra 
g that associates to any generator X a 6 g a parameter dependent generator 
X a (A) G Q (g) satisfying the following commutation relations 



[x* (A) , x? M] = £ rf x ^ x l~ XJ ^ . (3) 

7 

Here are the structure constants of the Lie algebra g; A and \i are complex 
spectral parameters. One can form Hermitian operators by using the dot 
product defined via the g-invariant metric tensor g a/3 ~ Tr (ad(X a )ad(X^)): 

K(X)=X(X)-X(X). (4) 

K (A) are not Casimir operators because they do not commute with all gen- 
erators of Q (g). However, these operators commute among themselves 

[K(\),K(n)]=0. (5) 

This implies that the system is integrable and K(X) are the integrals of 
motion. 

There is a unique rational realization^] of the generators X a (A) in terms of 
the generators of L copies of the algebra g given by the following expression: 

L 1 

r ( A ) = E7TT 1 " + ^ ( 6 ) 



i=i 



The Zi are L arbitrary numbers, which will ultimately be related to the 
single particle energies. Here we deviate from [16] by introducing the set 



2 There are non-rational Gaudin algebra realizations but we will not consider them here. 
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of arbitrary parameters p a . However, most of the expressions related to 
the Gaudin algebra (EJ) derived for the case p a = are still valid because 
X a (A) — > X a (A) + p a is an algebra isomorphism. The shift of the elements 
of Q by the p parameters is a key to the symmetry breaking in the model. 

3.1 Richardson-Gaudin Operators 

For the realization ([6]) of the Gaudin algebra, the integrals of motion are: 



C*2 is the second Casimir operator of the i-th. copy of g. Thus the first 
two terms are constants for states built on a tensor product of irreducible 
representations of q. The Ri are the Richardson-Gaudin operators [HI [2] and 
are one-half of the residue of K(\) at A = Z{ : 



By taking the residues of (jSJ at A = ^ and p = Zj one can see that the 
Ri operators commute among themselves. Therefore, they define a new set 
of integrals of motion. Thus any function of the Ri operators can be used 
as a model Hamiltonian for this integrable system. In particular, any linear 
combination of the Ri operators is at most quadratic in the generators. 

For a singular semi-simple algebras, the eigenvalues of the Ri operators 
can be obtained from the eigenvalues k(X) of the -^(A) operator given in [TB] 
by taking the appropriate residue (1/2 Res(fc(A), A = Zj)): 



Here e a , a are solutions of the generalized Richardson equations [161 131 EE] : 










5 



L is the number of copies of the algebra g and A, is the weight for the i-th 
copy. The Aj • ir a are actually the eigenvalues of the generators H a = n^h s 
at the lowest /highest weigh state of the i-th copy of g. The 7r" are the 
components of the positive simple roota^l ir a of the Lie algebra g in the Cartan- 
Weyl basis ([h s , E a ] = 7CgE a ), and £-7r a = 7r"p s , where p s are the components 
of the symmetry breaking one-body operator £ = p s h s along the Cartan 
generators h s . The rank of g is r and M a are positive numbers related to the 
eigenvalues m a of H a at the desired eigenstate (M a = Yli-^i ~ m i)- 

3.2 Symmetry Breaking 

Although one can use any function of Ri as a Hamiltonian, there is a par- 
ticular liner combination of the Gaudin operators Ri that results in a simple 
expression which is linear in the spectral parameters e a ^ a . From this expres- 
sion one can see that the breaking of the g-symmetry is due to the £ terms: 

h = j2 ZiRi = E *£i - ~ ( c a - c n ■ ( 10 ) 

i i \ i / 

In general zi ^ Zj, thus the term Yli z i£i would mix different irreps because 
it will not commute with the total Casimir operator C2 that is built from 
the generators X = Yli-^i- The final symmetry of this Hamiltonian is de- 
termined by the set of generators that commute with £. In order to see that 
the eigenvalues of (ITU]) are linear in the spectral parameters e ajQ one has to 
multiply the generalized Richardson equations (EJ) by e ajQ and sum over all 
the indexes. After some manipulations one would observe that the last term 
in OH]) appears in the relevant expressions. 

As we already discussed in the previous section, common feature of the 
dynamical symmetry nuclear models is that they are all defined for degener- 
ate single particle levels. Single particle energy splitting breaks the dynamical 
symmetry but may still preserve the integrability. The pairing model with 
non-degenerate single particle levels, whose exact solution has been found 
by Richardson, along with the above discussed Richardson- Gaudin construc- 
tions are examples of exactly solvable models with such characteristics [14]. 

Models based on fermion realization of the generators of the type a+a^", 
such as the sp(2n) algebras, are naturally suitable for non- degenerate sin- 

3 E a is a simple root vector if it cannot be witten as a commutator of any other two 
positive root vectors. 
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gle particle systems. This is related to the fact that the fermion realiza- 
tions of the corresponding Cartan generators are related to the fermion num- 
ber operators. This is an easy observation if one looks at the commutator 
[X^g,X~J = n a + np — 1 where the non-Cartan generators are in the form 
X*p = «„aj and (a a ) + = a+. Note that n a and np enter on an equal footing. 
In contrast models that are built on generators of the form X^g = a+ap, such 
as 5ll(n) algebras, result in Cartan generators of the form n a — rip. 

4 Pairing in Nuclei 

Now we will direct our discussion towards the applications of the general- 
ized Richardson-Gaudin models in nuclear physics. For this reason we shall 
specify Zj to be related to the single particle energies £, by the expression 
Zi = 2£j. p a to be non-zero only for X a that are elements of the Cartan 
algebra of g, that is p a ^ and a = l...r. Then a "pairing" like Hamiltonian 
Hp is obtained from (fTOj) by considering Hp = gH and its form is: 

Hp = j2 2£ * 5 - h i-9j2Y, Y ^ Y M 

r M a 

Ep = ^2^2e a<a S a , g=Y[—, 5 a = gp a - 

a=l a=l p a ^0 P 

The Ep are the eigenvalues of Hp and Yp are the positive root vectors with 
respect to the chosen Cartan algebra {h a } of g. Now g plays the row of a 
coupling constant for the two-body interaction term Y~?Yp. 

This is particularly clear in the case of the pairing model where g=su(2) 
with generators Y + = a^a^, Y- = a^a^ and h = a^a^ + ata^. Since there is 
only one Cartan generator h, there is only one p — ^, and 5=1, and Ep 

results in the usual expression of a sum over the pair energies E P = 52 i=1 e% 
for the standard pairing Hamiltonian ([T]). 

4.1 T=l Proton-Neutron Pairing as so (5) RG-model 

To better illustrate the current framework we will briefly discuss the T = 1 
pairing in proton- neutron systems which is related to an so (5) Lie algebra 
[13]. In this case the one-level system is constructed from various proton- 
proton, neutron-neutron, and proton-neutron pairs. By choosing the Cartan 
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generators to be the total particle number operator h 2 = 1 — (N p + N n )/2 
and the third projection of the isospin h± — T Q — (N p — N n )/2. We find 
the positive root vectors^ of the algebra to be the hard boson annihilation 
operators {&(//) : \i = 1,2,3} = {n^n^, (p^n^ + n^p^ ) /y/2, p^p± } plus the 
isospin rising operator T + = (p^n^ + p^n^) / ' \pl. The simple root vectors 
are {6(3), T + } and 6(3) = p~^p~l is the singular root vector [16]. In the chosen 
basis the corresponding commutation relations are: 

[h 2 M = 0, [&M,6(f)]=0, fM, v=l,2, 3 

[h 2 ,b(fi)} = b(ji), [h 2 ,T + } = 0, [h u T+]=T + 

[/n,6(l)] = 6(1), [/n, 6(2)1 = 0, [/n, 6(3)] = -6(3) 

[6(3),T + ] = 6(2), [T + , 6(2)] =6(1), [T + ,6(l)] = 0. 

By including the conjugated operators one closes the algebra 50 (5). The 
isospin sub-algebra sut(2) C so(5) is generated by the so(5) generators that 
are not related to the singular root vector 6(3) = pZpT- That is, su<r(2) is 
generated by T + , T_ = (T + ) + , and [T + ,T_] = T — h\. Even more, all the 
given so (5) generators can be recognized as sut{2) tensor operators. 

Since so (5) is a rank two algebra, we have two types of spectral param- 
eters: ei }a and e 2 ^ which we will denote by w a and vp in the following 
discussion. The upper bounds M 1 and M 2 , for the indices a and (3, are 
related to the isospin T and the total number of pairs M via the expressions 
M 1 = M — T and M 2 = M. The scalar products of the simple roots are 
(7r 2 ,7r 2 ) = 2, (7Ti,7Ti) = l,(7r 2 ,7Ti) = —1. If we consider now the spherical 
shell model, where protons and neutrons can occupy single particle states 
with quantum numbers (j,rrij) then the sub-index j" corresponds to rrij > 
and I to rrij < and the single particle index i labels the states (j, \rrij\). Due 
to the rotational symmetry, we can use the angular momentum j instead of % 
but have to take into account the corresponding degeneracy Qj = (2j + 1)/2. 
Finally the weights A° are the same for any % and correspond to the funda- 
mental representation of so(5), that is, A 1 = 0, and A 2 = 1. Putting all this 
together with the choice p 1 = 0,p 2 = —1/g, and Zi = 2Ei in the general- 
ized Richardson equating we obtain the equations for the proton-neutron 

4 Positiveness of a root vector is determined by the positiveness of the corresponding 
eigenvalues: first with respect to h r , if zero then one has to look at h r —i and so on. 
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T = 1 pairing that were given also by Links et al [T3] and Asorey et al [T] : 



- = Ett-^ + E^— + ( n ) 

g ^ 2ei- v a v a - Vp ^ w-y- v a 

i=l 0^=a K 7=1 ' 

M . M-T . M 

= y +y , £ = Vv 

a=l ' "7^7 a=l 

The spectral parameters v a have the same meaning as pair energies. If one 
allows for isospin breaking then one has to set p 1 = A ^ [3]. 

4.2 Spin-Isospin pn-pairing as so (8) RG-model 

In the previous section on our discussion of the T=l proton-neutron pairing 
as so (5) RG-model, we considered the intrinsic symmetry space to be the 
isospin SUt{2) symmetry while the extrinsic spaces labeled by % were related 
to the total spin states (j, rrij). In particular, the pairs operators {b(fi), b + (fi)} 
and the sut(2) algebra generators {T Q ,T±} were time-reversal invariant. For 
the study of the spin-isospin pn-pairing it is appropriate to consider the 
Wigner's SUst(4:) = SUs(2) x SUt(2) as intrinsic symmetry of the system 
[TT] and the orbital angular momentum (I, mi) as the extrinsic space labeled 
with i in the appropriate sums. This way the protons and neutrons are 
described by the operators: ai ijms>a .^ T where s — t — | or briefly ai um ^ T . 
The pair operators are defined as isovector and spinvector tensor operators: 

Pri = Ke^JooV ' D <n = V^h [ a h a k]ol°o ■ 

The 1\ = (2/j + l)/2) appears here due to the structure of the Clebsch- 
Gordan coefficients < lm,l — m\00 >. The Wigner's su(4) along with the 
u(l) number operator N = N nt +N n ,+N Pt +N p , = 2n which is double of the 
pair number operator n, form w(4) algebra with the following 16 generators: 



X Tl0 - lT2 a 2 i ^ ] 1 



rn,Ti<Ti a km,T2<J2y 



that are u(l), 5115(2), and sut(2) tensors up to a factor yH^: 

at r + i 000 c r + 1 010 t r + 1 001 v r + 1 011 
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The relevant hamiltonian has equal spin and isospin pairing strength: 

L 



13,11 



with eigenvalues: E = J^^ 1 e « where each Mi pair contributes pair energies 
e a determined by the four equations [T2] : 

, L n Ah M 2 1 

1 Sij z 



Mi 1 M 2 „ M 3 M 4 



a' a'(^a) a' or 

M 2 1 M 3 











1 




10a! 




Va 




1 





U) a i 
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Va' 








2 




la' 




7a' 



a'(7=a) 
M 2 M 4 

a' a('^=a) 

Since this is 50(8) RG model of rank is 4 there are 4 sets of spectral pa- 
rameters. The number of spectral parameters in each set is determined 
by the relevant u(4) Wigner multiplet. For a given number of pairs M, 
these w(4) multiplets can be classified using Young tableaux. Each multi- 
plet is defined by a partition of M in 4 numbers, [A1A2A3A4], constrained by: 
Yli^i ^ Ai > A2 > A3 > A4 > 0. The labels Aj are related to the number of 
pairs in the lowest /highest weight state. Thus Mi = Ai + A2 + A3 + A4 is the 
number of pairs M, M 2 = A2 + A3 + A4, M 3 = A3 + A4, and M 4 = A4. 



5 Conclusions and Discussions 

Using the Cartan classification of the semi-simple Lie algebras one can see 
that many physics models can be generalized within the above framework. 
In Table [5] are shown the semi-simple Lie groups up to rank four and the 
corresponding main fermion models. R&R denotes the model developed by 
G. Rosensteel and D. J. Rowe and TCF stands for Trapped Cold Fermions 
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rank n 


A n su(n + 1) 


B„ so(2n + l) 


C„ sp(2n) 


D n so(2n) 


1 


su(2) pairing 


so(3) ~su(2) 


sp(2) ~i3u(2) 


so(2)~u(l) 


2 


su(3) Elliott 


so(5) T=l pairing 


Sp(4) ~ so(5) 


so(4) ~ su(2) ©ssu(2) 


3 


su(4) Wigner 


so (7) C so(8) FDSM 


Sp(6) R&R 


so(6) ~ su(4) TCF 


4 


su(5) 


50 (9) 


sp(8) 


So(8) Evans, FDSM 



Table 1: Group structure associated with important nuclear physics models; 
pairing models are in boldface. Isomorphisms of Lie algebras are denoted 
with ~. 



The RG-models are a new important mathematical tool to study the 
behavior of physical systems. They are finding applications to the such fields 
of studies as super- conducting grains, atomic nuclei, and trapped fermion 
atoms. 
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